m1=1; v1=[0 1];
m2=1; v2=[8 1];
m3=1; v3=[2 4];

m=m1+m2+m3;
v=1/m*(m1*v1+m2*v2+m3*v3)

A=[1 1 1; 0 8 2; 1 1 4]
b=[6 8 12]'
w=A\b
A*w

m=m+6;
v=1/m*((m1+w(1))*v1+(m2+w(2))*v2+(m3+w(3))*v3)
